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ABSTRACT 

The  numerical  frontogenesis  model  of  Williams  (1973) 
is  modified  to  include  moisture  with  its  subsequent  con- 
densation and  release  of  latent  heating.   The  turbulent 
diffusions  of  momentum,  heat  and  moisture  are  repre- 
sented with  various  coefficients.   The  numerical  solutions 
show  realistic  quasi-steady  fronts  forming  within  one  to 
two  days.   These  solutions  are  examined  and  compared  over 
a  range  of  the  various  coefficients,  and  various  ranges  of 
temperature . 

Inclusion  of  moisture  in  the  model  causes  intensifi- 
cation of  baroclinicity  at  mid  and  upper  levels.   Also 
noted  is  a  sensitivity  of  the  moist  model  to  reference 
potential  temperature  due  to  the  exponential  relationship 
between  saturation  vapor  pressure  and  temperature. 
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I.   INTRODUCTION 


Hoskins  and  Bretherton  (1972)  have  shown  analytically, 
and  Williams  (1967,  1972)  has  shown  numerically  that  dis- 
continuous fronts  can  form  within  a  finite  period  of  time 
if  no  turbulent  -diffusion  is  present.   These  studies  sug- 
gest that  a  discontinuity  will  form  within  24  to  36  hours 
when  reasonable  initial  conditions  are  used.   If  turbulent 
diffusion  is  present,  it  can  be  expected  that  a  balance 
will  be  achieved  between  the  f rontogenetic  advections  and 
turbulent  diffusions  of  momentum,  heat,  and  moisture.  The 
front  should  remain  in  this  state  of  quasi-balance  as  long 
as  the  large  scale  deformation  field  causes  f rontogenetic 
advections  around  the  front. 

The  numerical  f rontogenesis  model  of  Williams  (1973) 
shows  the  effects  of  horizontal  and  vertical  turbulent 
diffusions  of  momentum  and  heat  in  the  formation  of  these 
quasi-steady  fronts.   This  model,  which  also  includes  an 
Ekman  boundary  layer,  will  hereafter  be  referred  to  as 
W73  and  the  model  as  it  existed  in  Williams  (1972)  will  be 
referred  to  as  W72. 

This  investigation  appends  moisture  with  its  subse- 
quent condensation  and  release  of  latent  heating  to  W73. 
It  is  assumed  that  all  the  latent  heat  released  in  this 
condensation  is  used  to  heat  the  air  isobar ical ly . 


The  purpose  of  this  study  is  to  obtain  and  examine 
frontal  solutions  from  the  model  throughout  a  four  day 
period.   The  physical  model  is  essentially  the  same  as 
one  of  the  models  treated  by  Hoskins  (1971)  and  Hoskins 
and  Bretherton  (1972).   The  basic  model,  as  it  now  exists, 
numerically  describes  f rontogenesis  which  is  forced  by  a 
nondivergent  horizontal  wind  field  which  contains  stretch- 
ing deformation.   This  deformation  wind  field  is  constant 
in  time  and  is  independent  of  height,  except  in  the  Ekman 
boundary  layer.   The  moist  hydrostatic  primitive  equations 
with  diffusion  are  used  in  which  the  time  dependent  quan- 
tities are  functions  of  Y  and  P  only.   Therefore,  output 
from  the  model  is  in  the  form  of  vertical  cross-sections 
in  the  Y-P  plane. 

W73  has  been  transformed  to  pressure  coordinates  to 
obtain  more  accurate  values  of  pressure  for  the  moisture 
equations  that  have  been  included  in  the  present  model. 

In  Section  2,  the  basic  forecast  equations  are  devel- 
oped,  transformed  to  pressure  coordinates  and  a  simplifi- 
cation which  keeps  the  problem  two  dimensional  is  discussed 
The  initial  conditions  are  presented  and  transformed  to 
pressure  coordinates  in  Section  3.   Convective  adjustment 
processes  and  latent  heat  release  in  the  model  are  included 
in  Section  4.   The  experiments  undertaken  are  examined  in 
Section  5  and  conclusions  are  given  in  Section  6. 


II.   BASIC  EQUATIONS  AND  BOUNDARY  CONDITIONS 

The  hydrostatic  primitive  equations  with  diffusion 
and  moisture  may  be  written  as  follows: 

3V  3(wV)  32V 

-^  +V'(W)+   -     +  V<J)  +  f  (KxV)=A  V2V+C ^    ,  (2.1) 

3s 

||  +v.(ev)+  ^(we)-  AQv2e+ce  i!|  +  qa+  qm  ♦  hc  ,     (2.2) 

3s 

I?  tv-(«y)+  |j(w«)-  V2q+Cq  ff  +  ma  "  Mc  '  (2'3) 

3s 

V'V   +  ~  =  0  ,  (2.4) 

~     ds 

■P-  =  L(s)9  ,  (2.5) 

ds 


where  the  vertical  coordinate  has  been  replaced  by  s  to 
aid  in  later  discussions  on  transformation.   The  equations 
are  from  W73  with  slight  modification  due  to  the  addition 
of  the  equation  of  conservation  of  water  wapor  (2.3)  which 
is  analogous  to  the  thermodynamic  equation  (2.2).   The 
function  L(s)  in  the  hydrostatic  equation  (2.5)  is  a  con- 
stant g/8   in  height  coordinates.   Also  sinks  and  sources 
o 

of  heating  and  moisture  have  been  added  to  the  thermody- 
namic equation  (2.2)  and  the  conservation  of  water  vapor 
equation  (2.3)  . 

The  notation  used  is  basically  the  same  as  W72  which 
uses  height  as  the  vertical  coordinate.   However,  in  our 
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subsequent  transformation  to  pressure  coordinates,  0  is 

defined  as  a  departure  from  the  domain  average  of  the 

initial  state,  <f>  represents  geopotential ,  gz  and  w  is 

ds/dt  in  the  continuity  equation  (2.4),  thereby  allowing 

compressibility. 

The  quantities,  A  ,  A0,  and  A   are  the  horizontal  tur- 

mo        q 

bulent  diffusion  coefficients  of  momentum,  heat,  and 

moisture,  respectively,  while  C  ,  CQ,  and  C   are  the  cor- 

e  m    6        q 

responding  vertical  coefficients.   These  constant  coeffi- 
cients were  assigned  various  reasonable  values  in  the 
different  moisture  experiments.   Dry  and  moist  convective 
adjustment  processes  are  represented  by  the  functions  Q 
and  Q   while  latent  heating  due  to  large  scale  condensa- 
tion is  represented  by  H  .   In  the  moisture  equation,  M 
is  a  sink  or  source  of  moisture  due  to  moist  convective 

adjustment  while  M„  is  a  sink  for  moisture  due  to  large 
J  C 

scale  condensation. 

Boundary  conditions  for  the  above  equations  are  such 
that  the  domain  is  bounded  by  two  rigid  quasi-horizontal 
planes.   The  no-slip  condition  at  the  lower  boundary  and 
the  no-stress  condition  at  the  top  roughly  represents  the 
troposphere.   The  conditions  are  as  follows: 


w(x,y,o,t)  =  w(x,y,H,t) 
V(x,y,o,t)  =  0 

Cm  37<**y'H,t) 

ce  TF(x'y'°'t) 

c„  ^7-(x,y,o,t) 
q  os 


=  0 


=  0 


,   (2.6) 


ae 


Ce  _(x,y,H,t) 

cq  if <*.y.H,t) 


=  0 


=  0 
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where  H  represents  the  distance  between  the  upper  and 
lower  boundaries.   The  lower  boundary  conditions  are  realis- 
tic for  fronts  over  level  terrain.   Therefore,  fronto- 
genesis  in  the  lower  troposphere  is  more  realistic  in  the' 
model  than  are  conditions  near  the  " tropopause" .   The  addi- 
tional boundary  conditions  for  moisture  are  completely 
analogous  to  the  heating  conditions  in  that  fluxes  of  the 
quantities  do  not  pass  through  the  boundaries. 

The  approximate  steady  state  solutions  to  the  hydro- 
static primitive  equations  from  [  (2 . 1) -  (2 . 5) ]  are: 


V  =  U  =  D{[x(l-e    cos  BS)+  ye    sin  Bsl  i  ,\ 
•v    «*  *» 

-Be  -Bs  i 

-  [y(l-e    cos  Bs ) -  xe    sin  Bs]  j ) 

w  =  0  , 

2   2     2 

<(>  =  $  =  -D  (x   +  y  )/2  -  fD  xy  +  F(s) 

6  =  G(s) 
q  =  H(s) 


(2.7) 


1/2 
where  B  =  (f/2C  )     ,  D  and  f  are  constant  and  where  F(s) 

m 

is  such  that  the  hydrostatic  equation  (2.5)  is  satisfied. 
The  functions  G(s)  and  H(s)  are  arbitrary  functions  of  s 
since  w  =  0.   These  solutions  satisfy  the  boundary  condi- 
tions (2.6)  and  allow  for  inclusion  of  the  Ekman  boundary 
layer.   These  relations  (2.7)  satisfy  the  steady-state  equa- 
tions if  some  small  advection  terms  in  the  boundary  layer 
are  neglected.   These  small  terms  are  usually  dropped  in 
the  Ekman  Theory. 
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The  equations  and  the  model  become  two  dimensional 
after  a  subdivision  of  the  dependent  variables  as  follows 


V  = 


w  = 


e 
q 


U(x,y,x)  +  u(y,s,t)i  +  v(y,s,t)j 

*v  ****** 

w(y,s, t) 

6(y,s,t) 

$(x,y)  +  IT  (y,s,t) 

q(y,s,t) 


(2.8) 


It  is  noted  that  all  departures  from  the  approximate 
steady  deformation  solutions  (2.7)  are  assumed  to  be  in- 
dependent of  x.  If  we  substitute  the  relations  (2.8)  in 
the  i  component  of  (2.1)  we  obtain: 


dU  9  /  X  3  /  N  rWI  "BS  T,         X 

■srr   +    i\—    (uv)     +   t-(wu)     +    uD(l-e  cos    Bs) 

dt  dy  ds 

~,  ,,    "Bs     „  .     -Bs     „  ,3u     „  -Bs  . 
-D[y(l-e    cos  Bs)-xe    sin  Bs]t —  +  vDe    sin  Bs 

3y 


3      -Bs  —  B  s 

+  wD  ■= — [xe    (1-cos  Bs)+ye    sin  Bs]=fv  +  A 
d  s  m 


92u 

3y2 


+  c 


a2 


u 


m  3s2 


(2.9) 


Note  that  the  assumption  that  u  is  independent  of  x 
for  all  time  is  violated  in  (2.9)  in  the  two  terms  where 
x  appears.   This  indicates  that  the  other  quantities 
which  were  assumed  independent  of  x  are,  in  fact,  x- 
dependent  and  this  would  bring  other  terms  into  the  equa- 
tions.  In  this  study  we  will  neglect  these  added  terms 
and  apply  (2.9)  and  the  other  equations  at  x  =  0.   If  the 
initial  ufv#w,9,Tr  and  q  fields  are  independent  of  x,  the 
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error  from  this  approximation  will  grow  slowly  and  will 
be  confined  to  the  boundary  layer  since  the  x-dependent 
terms  in  (2.9)  are  zero  outside  the  boundary  layer. 
Since  f rontogenesis  occurs  very  rapidly  it  is  expected 
that  the  development  of  x-variations  in  the  dependent 
variables  would  have  only  a  small  effect  on  the  resulting 
quasi-steady  front.   In  any  case  these  effects  could  not 
be  observed  in  the  atmosphere,  because  atmospheric  fronts 
have  some  variation  in  the  basic  fields  along  the  front 
which  would  be  much  more  important. 

When  equation  (2.9)  is  evaluated  at  x  =  0  it  becomes 


3u    3(uv)   3(wu)  r  3v  .r„3u    v3v,  -Bs  . 

■arr  +  -k—    +  -5—    -Iu^—  +TV-5—  -  — 5— (e    sin  Bs  ) 

3t    3y       3s        dy      dy      dy 


3  ,  -Bs  .     .        32u 

-Vw  tt— (e    sin  Bs)=fv+A   +  C 

ds  m  „  2     m 

3y 


32u 
3s2 


(2.10) 


Bs 
where  V  =  -Dy  and  T    =  1-e   cos  Bs .   The  j  component  of 

(2.1)  applied  at  x  =  0  takes  the  form 


-~-r   +   -5— (vv)+  t-  wv  +Tv-vV  -uv-  e    sin  Bs 
3t    3y       3s       3y       3y 


-(wV)|-(e"B  cos  B  )=-|^  -fu+A 
ds  dy       m 


32u 


3y 


+  C 
2    m 


32v 
3s2 


(2.11) 


In  a  similar  manner  from  equation  (2.2),  the  following 

equation  for  the  departure  of  potential  temperature  can  be 

derived : 

36  ,  3(v6)    3(w6)    p„  39    B   320  ,  „   329 

■grr    +    — 5 +  —5 +  TV  -5—  =  A   +  CQ  

3t     3y       3s         3y     8  ^..2     9  ^_2 


3y- 


3s' 


+  Q   +  Q   +  H 
a     m     c 


(2.12) 
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and  the  following  equation  for  q  follows  from  (2.3), 


3q    3(vq)     3 


3q   _ 


3t     3y 


+  ^-(wq)  +  TV  -^ 
3s   ^        3y 


a   i!a  +  c   i!a 


q  3y2 


3s- 


+  MA  -Mc  , 


Equations  (2.4)  and  (2.5)  become 


(2.13) 


3v    3w 
3y    3s 


"577  +  ^r  -  0  , 


and 


3tt 

3s 


=  l(s)   e  . 


(2.14) 


(2.15) 


The  boundary  conditions  at  the  bottom  and  top  of  the  domain 
respectively  are 


u  =  v  =  w  =  C 


36 


=  C 


3q 


6  3s     q  3  s 


3u 


'm  3 1 


„   3v 

C    TT—  =  w 

m  ds 


36 


0  ,  s 

r,     3q 


=  0 


'6  3s     q  3s 


s  =  H 


(2.16) 


If  we  define  the  vertical  average  of  a  quantity  as 

H 


<<)>-£/   (   )  ds  , 
n 

o 


(2.17) 


integrate  the  hydrostatic  equation  (2.15)  with  respect  to 
s,  and  remove  the  vertical  mean,  it  follows  that 


s  s 

T7-<ir>  =  [  /   L(s)8ds-  <  /  L(s)6ds>] 
o  o 


(2.18) 


Take  the  vertical  average  of  (2.14)  and  use  the  boundary 
(2.16)  which  gives 

3 


3y 


<v>  =  0  . 


(2.19) 
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This  equation  states  that  the  total  mass  flux  in  the  j 
direction  is  independent  of  y.   In  W72  a  symmetry  argument 
was  used  to  show  that  this  flux  must  be  zero.   This  argu- 
ment does  not  hold  strictly  in  the  present  experiment, 
but  it  can  be  expected  that  at  large  distances  from  the 
frontal  zone  the  disturbance  mass  flux  will  vanish.   Thus, 
it  is  assumed  that 


<  v  >  =  o  . 


(2.20) 


If  the  development  of  W72  is  used,  equation  (2.11)  can 
be  written 


f^  +  |-[vv  -<vv>]+  |~(wv)  +  |~(TvV  -  <Tv>  V) 
dt    dy  ds         dy 

dV.   -Bs     „   ^   -Bs     „  v   „,   8   .  -Bs     „  , 
-  7; — [ue    sxn  Bs-<ue    sin  Bs>  -V[w  7; — (e    cos  Bs) 
dy  ds 


-  <w  v~(e   Scos  BS)>]  =  -  -s—  (IT  -<ir>) 
ds  dy 


-f (u-<u>) 


+  A 


m 


ifv  +  c  fi!v  +  1  |v 


dy 


"     ds2 


H  3s  s=0 


] 


(2.21) 


The  equations  (2.10),  (2.12),  (2.13),  (2.14),  (2.18) 
and  (2.21)  form  a  complete  set  which  can  be  solved  by  a 
pure  marching  process.   The  finite  difference  equations 
conserve  mean  squares  in  the  advection  terms  and  are  de- 
scribed by  Williams  (1967).   The  Matsuno  scheme  is  used 
everywhere  except  in  the  arrangement  of  variables  and  the 
finite  difference  approximations  are  the  same  as  those 
used  by  Williams  (1967).   In  order  to  close  the  problem, 
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computational  boundaries  must  be  introduced  in  y.   Since 
the  disturbance  velocities  should  die  out  at  a  sufficient 
distance  from  the  axis  of  dilation,  then 


v(+  Y,s,t)  =  0  . 


(2.22) 


However,  there  is  appreciable  inflow  across  these  com- 
putational boundaries  since  V(+^  Y)  =  +_   DY  .   The  quantities 
uf  0,  and  q  which  are  advected  across  the  boundaries  must 
be  specified  independent  of  the  interior  values  if  compu- 
tational stability  is  to  be  maintained  (Platzman,  1954). 
Thus,  the  following  boundary  conditions  are  used: 


u[+(Y  +  Ay/2),s,t]  =  u  [+(Y  +  Ay/2),s,0] 
9[+(Y  +  Ay/2);S,t]  =  6  [+(Y  +  Ay/2),s,0] 
qt+(Y    +    Ay/2),s,t]     =    q     [+(Y    +    Ay/2),s,0] 


(2.23) 


The  computational  boundaries  y  =  +  Y  are  placed  between 
gridpoints  so  that  the  above  conditions  are  actually 
applied  at  y  p  +  (Y  +  Ay/2). 

The  boundary  conditions  were  found  to  be  satisfactory 
except  that  a  solution  separation  developed  near  the 
boundaries.   This  was  controlled  by  using  a  Matsuno  time 
differencing  scheme  for  stability. 

In  order  to  have  precise  values  of  pressure  at  grid- 
points  for  the  moisture  experiments  it  was  deemed  neces- 
sary to  transform  W72  to  pressure  coordinates. 
Minimization  of  changes  to  finite  differencing  schemes 
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of  the  model  have  been  achieved  by  defining  a  parameter  s 
which  is  analogous  to  a  field  of  z.   Specifically, 

s  =  Pq  -  P  ,  (2.24) 

a  field  that  will  increase  in  the  upward  direction  and 
from  which  values  of  pressure  can  easily  be  obtained. 
Further,  it  is  seen  that  9s  =  -  3P  so  that  finite  differ- 
encing schemes  in  W72  can  be  equally  well  applied  to  the 
field  s. 

Through  the  expression  in  s  coordinates  the  momentum 
equations  (2.10)  and  (2.21),  the  thermodynamic  equation 
(2.12)  and  the  moisture  equation  (2.13)  remain  as  before 
for  the  model.   However,  the  continuity  equation  (2.14) 
in  the  model  actually  remains  in  w  and  describes  a  nega- 
tive field  of  omega  in  that, 

-  ds  ,  _  _  _ , 

w  =  — -  =  -0)  ,  (2.25) 

at 

Equation  (2.5),  the  hydrostatic  equation,  differs 
significantly  due  to  the  transformation.  The  derivation 
is  briefly  as  follows.   The  hydrostatic  equation  expressed 
in  pressure  coordinates  and  the  equation  of  state  gives 

|i  -  -O  -  -  2E  .  (2.26) 

9p  p 

If  a  form  of  the  potential  temperature  definition  is 

used  , 

T  =  6 (£-)K  ,  (2.27) 

Po 
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and  combined  with  (2.26)  and  (2.27)  it  follows 


3P      P   P, 


o 


It  can  be  shown  that 


S  <?->K=  C   ^(^)K   .  (2.29) 

P   PQ      P  9P  PQ 


The  latter  form  simplifies  finite  difference  energy  conser- 
vation.  Therefore,  substituting  equation  (2.29)  into 
equation  (2.28)  and  transforming  to  s  coordinates  we  obtain 


ds        p  ds    p 

r  o 


a  form  of  the  hydrostatic  equation.   Note  that  L(s)  is  not 
constant  as  in  height  coordinates  since 


L(s)  =  —(£—)       =  -  C   -5—  t^1— )   .  (2.31) 

p  Po         p  3s  pQ 


It  follows  that  the  complete  set  of  equations  solved  by 

the  marching  process  remains  intact  but  with  L(s)  defined 
as  in  (2.31) . 

The  boundary  conditions  (2.16)  also  remain  intact 
except  that  the  lower  boundary  of  vertical  motion  in  the 
model  is  approximated  as 

w=-w=0,ats=0.  (2.32) 
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III.   INITIAL  CONDITIONS 

The  initial  conditions  in  height  of  s  coordinates  are 
the  same  as  those  employed  in  W73.  The  initial  potential 
temperature  field  is  given  by 


89I      H       2 
6(y,z,o)  =  [-r (z-  — )  ]  -a  (— )  arctan  (sinh  a  y) 

where 

a  =  .  (3.1) 


4 


36i 


The  quantity  38  /9z  which  is  constant  is  the  initial  static 
stability,  8   being  the  initial  horizontally  averaged 
initial  theta  and  "a"  is  one-half  the  total  horizontal  tem- 
perature variation. 

By  substituting  equation  (3.1)  into  the  thermal  wind 
equation  and  integrating  we  obtain  the  x-component  of  the 
rotational  wind  which  is 

,       ,     2  gaa  .      H.  ._  .  . 

u(y,z,o)  =  -  2_—  (z  -  — )  sech  a  y  .  (3.2) 

U  ID  2 

o 

The  divergent  part  of  the  wind  or  v  field  is  obtained 
from  the  quasi-geostrophic  equations  in  z  coordinates. 
Details  are  as  in  Williams  (1972).   The  remaining  fields 
of  w  and  it  are  obtained  from  the  continuity  (2.14)  and 
hydrostatic  (2.15)  equations. 


20 


Once  again  transformation  of  the  above  initial  fields 
must  be  made  to  pressure  coordinates. 

To  obtain  initial  conditions  for  the  theta  field  we 

must  first  find   0   ,  the  horizontally  averaged  initial 

theta  based  on   0,.  (p  )  .   Derivation  is  as  follows.   Given 

I   o 

knowledge  of  differentials  and  the  definition  of  geopoten- 
tial  it  is  seen  that 

30T   30T  a  .    aeT    .. 

3p     3z   3P  "  g  3z   3p  '  *J,J; 

where  30  /3z  is  the  prescribed  static  stability  at  t  =  0. 

Substitute  a  definition  of  potential  temperature  (2.27) 
to  obtain 


rr-Z-   =  -5-i  t-  -  -  (E-)*]  0T  .  (3.4) 

dp     dz      g  p   p       I 


Integrate  from  p   to  some  p  , 

rpi  35i  „  rp*  ^i  i.p  .k  a 

P„      I  P^ 

O  O 

which  yields 

0              30  p 

1„  JL              --**-*  EC*-)*  "  (-^)K1   •               (3.6) 

0l(Po)       g  8Z   K   Po  Po 


where   K  =  R/C   .  it  follows  that 
P 

C   30 

0T(P)  =  0T(P  )  EXP [-£  *-±-(l- C2-)  )]  .  (3.7) 

I         I   o       g   dz      p 

o 
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Equation  (3.7)  is  representative  of  term  one  of  equa- 
tion (3.1).  Therefore,  the  initial  potential  temperature 
field  in  pressure  coordinates  is  given  by 

C   39, 


O^y^o)  =  G^pJ  EXP[^^ 


i(l-(2-)K)] 
po 


-a(2/ir)  arctan  (sinh  ay). 


(3.8) 


where 
a  ~  — 


fTT 


f-rr 


4 


30. 


Hl/L(s  )  (-_— -) 
*    m    3p 


H 


1  p  p 
»  m   o 


)K(- 


3p   m 


(3.9) 


In  the  model  a  is  given  a  constant  median  value  by 
solving  equation  (3.10)  at  p  =  p   =  60  cb  and  using  the 
static  stability  30  /3p  at  this  level.   Recall  also  that 
pressure  in  equation  (3.8)  can  easily  be  obtained  for 
model  purposes  as  p  =  p   -  s. 

Further  it  should  be  recalled  that  9  is  a  departure 
from  the  domain  average  of  0  .   This  domain  average  is 
found  in  the;  model  by  finding  the  vertical  average  of 
equation  (3.7),  <  ©>• 

The  initial  x-component  of  the  velocity  or  initial 
u  equation  is  found  using  a  form  of  the  hydrostatic  equa- 
tion (2.26)  and  the  definition  of  geostrophic  wind, 


u  =  - 


1   3<|> 


(3.10) 


f   3y   ' 
A  form  of  the  thermal  wind  equation  is  found  to  be 


3u   _        13      3<j>   _    t    1    R,p_>k    3_6 

3p  f    3y    3p  f    p   p  3y 


(3.11) 
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If  we  again  assume  a  median  pressure  level  p   =60  cb, 

and  that  at  this  level  u(p  )  =  0,  we  integrate 

m 


P  P 

f  pL   dp  =  *  |i  f   I(p_)< 
J  3p     f  3y   /  p  p0 


dp 


(3.12) 


m 


m 


Recalling  that  u (p  )  =  0  and  that  K  =  R/C   the  integra- 

3         ^m  p  3 

tion  yields 

,      .     p  36.. p. K  ,pm.K..                   , ,  ... 

u(y,p,o)  =7^3—  [  (r— )  -  (— )  1  •                 (3.13) 

r   oy    p  p 

o  o 


From  (Williams,  1972)  equation  (4.3),  it  is  seen  that 


I6- 

3y 


2aa 


sech  (ay) 


(3.14) 


Therefore,  the  initial  u  field  in  pressure  coordinates  is 

given  by 

C  p 

P   2aa  .  ,p  ,K     .  m.K,     ,  ,  „  ,  ,. , 

u(y,p,o)  =  -  j^   — —  [  (f— )   -  (— )  ]  sech  (ay)  .    (3.15) 

o        o 

The  initial  divergent  wind  v  is  determined  as  in  Wil- 
liams (1972)  but  using  the  quasi-geostrophic  equation  in 
pressure  coordinates.   The  stream  function  is  defined  as 


and 


v  =  W/dp 
0)  =-d\p/dy 


(3.16) 


and  equation  (4.10)  of  Williams  (1972)  takes  the  form 

,2     m   dy 


9p' 


m 


3p   P 


(3.17) 


m  8y ' 
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where  L(s)  is  assumed 


L(s)  =[£  <^)K] 
Pm  Po 


The  equation  (3.17)  is  solved  for  ij>  from  the  initial 
temperature  field  (3.8)  with  the  technique  of  Ogura  and 
Charney  (1962)  and  v  is  then  determined  from  (3.16). 

It  is  noted  that  these  initial  conditions  do  not 
satisfy  the  boundary  conditions  on  u,  v,  and  0  at  the 
upper  and  lower  boundaries.   A  period  of  adjustment  will 
be  required  to  form  the  surface  friction  layer  and  some 
oscillations  may  be  observed  later. 

The  fields  of  w  and  IT  are  obtained  through  the  con- 
tinuity and  hydrostatic  equations  (2.14)  and  (2.15). 

Finally  it  is  necessary  to  define  explicitly  methods 
used  in  determining  the  initial  field  of  specific  humidity, 
given  Q  (y,p,o)  from  equation  (3.8). 

Saturation  vapor  pressure  e  ,  is  first  determined  at 
each  gridpoint  through  use  of  an  integrated  form  of  the 
Clausius-Clapeyron  equation,  an  equation  expressing  satura- 
tion vapor  pressure  as  a  function  of  temperature.   Accord- 
ing to  Hess  (1959) 

L 


e   =  .611  EXP  [-£■(£- 
s  R   T 

v   o 


,-n   , 


(3.18) 


where  e   is  expressed  in  centibars,  L   is  the  latent  heat 
s  c 

of  condensation,  R    is  the  gas  constant  for  water  vapor, 

and  T   is  a  constant  of  integration. 
o 
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Knowing  saturation  vapor  pressure,  the  saturated  speci- 
fic humidity  was  obtained  from  the  approximated  relation 
from  Haltiner  and  Martin  (1957) 


.622e 


qs  = 


(3.19) 


Therefore,  by  having  knowledge  of  potential  temperature 
at  t  =  0  from  equation  (3.8),  q   can  be  found  at  each  point 
using  equations  (3.18)  and  (3.19).   By  defining  q  to  be  a 
given  fraction  of  q   in  different  experiments,  atmospheres 
of  different  constant  initial  relative  humidities  were 
tested. 
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IV.   CONVECTIVE  ADJUSTMENT,  LATENT 
HEATING  AND  PRECIPITATION 


Inclusion  of  two  additional  heating  terms  in  the  moist 
thermodynamic  equation  (2.3)  necessitated  some  further 
amendment  to  W73.   Atmospheric  heating  by  dry  convective 
adjustment,  Q  ,  was  included  in  W73  to  eliminate  any  in- 
tense grid  scale  convection  which  developed  due  to  the  for- 
mation of  any  dry  unstable  lapse  rates.   This  redistribution 
of  heating  was  accomplished  using  averaging  techniques  to 
ensure  that  the  dry  adiabatic  lapse  rate  was  not  exceeded 

dp  — 

Lapse  rates  were  tested  in  each  vertical  column  of  the 

3  8 
model.   Upon  encountering  an  unstable  layer  (^ —  >  0),  a 

dp 

cumulative  average  was  taken  and  applied  to  the  gridpoints 

<-)  9 
within  this  layer  until  neutral  stability  (-~ —  =  0)  was 

dp 

attained.   Simultaneously  the  model  was  checked  for  insta- 
bilities developing  below  the  convectively  adjusted  layer 
and  potential  temperatures  at  these  levels  were  also  cumu- 
latively averaged  and  reassigned  the  new  values.   The  final 
result  of  the  technique  is  that  the  heating  is  redistri- 
buted to  simulate  physical  processes  of  the  atmosphere  and 

3  fi 
mathematically  that  7; —  <  0  throughout  each  column  of  the 

dp  — 

model.   A  completely  analogous  scheme  was  used  for  saturated 
parcels  in  the  moist  experiment,  only  using  equivalent 
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potential  temperatures,  0   ,  in  place  of  potential  tem- 
perature,  8. 

With  the  addition  of  moisture, consideration  of  large 
scale  condensation  and  the  subsequent  release  of  latent 
heating  became  important. 

According  to  Haltiner  (     ) ,  changes  in  moisture  and 
temperature  due  to  large  scale  condensation  were  deter- 
mined as  follows.   Where  relative  humidity  exceeds  100%, 
the  adjustments  in  temperature  6t  and  specific  humidity 
Sq  are  determined  by  the  following  equations: 


and 


q  +  6q  =  qg (T  +  6T,  p)  , 


C  6T  =  -  L   6q  . 
P  c 


(4.1) 


(4.2) 


These  equations  require  that  the  final  state  of  the 
air  be  exactly  saturated  and  that  excess  moisture  be  con- 
densed isobar ically ,  releasing  the  latent  heat  to  the  air 

Equations  (4.1)  and  (4.2)  can  be  expressed  in  itera- 
tive form  as 


3qs 

Aq  =  qs(l-r)  +  <^-)p  At  ;   q  =  rqs  , 


and        C  AT  =  -  LAq  . 
P 


(4.3) 
(4.4) 


Eliminating  AT  from  (4.3)  and  (404)  leads  to 
Aq  = 


qs   (1-r) 


L   8qs 

1  +  -  ( - 

c  13t 
p 


(4.5) 
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At  constant  pressure  the  Clausius  Clapeyron  equation 
can  be  expressed  as 


L   q 
s  %      c   s 


3q^ 


(4.6) 


Substituting  (4.6)  into  (4.5)  we  obtain 


Aq  = 


qs(l-r) 


A. 


(4.7) 


1  + 


CRT 
p  V 


Now  Aq  may  be  calculated  from  (4.7)  and  the  correspond- 
ing value  of  At  from  (4.4). 

The  first  estimates  of  the  adjusted  temperature  and 

specific  humidity  are,  therefore,  ■ 


T'  =  T  +  AT  and  q *  =  q  +  Aq  . 


(4.8) 


There  values  in  turn  were  used  to  obtain  improved  ap- 
proximations.  Iterations  were  containued  until  the  values 
of  T"  and  q*  at  each  level  matched  with  sufficient  accu- 
racy one  hundred  percent  relative  humidity.   The  changes 
were  then  applied  to  the  potential  temperature  and  speci- 
fic humidity  fields  of  the  model  accounting  for  the  term 

H   of  equation  (2.2)  and  M   of  equation  (2.3) 
c  c 

The  model  was  then  investigated  for  unstable  moist 
convection  and  its  associated  heat  of  condensation.   A 
convective  adjustment  scheme  similar  to  the  dry  one  was 
used  on  the  saturated  parcels  so  that  moist  static  stability 
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was  always  maintained.   A  field  of  equivalent  potential 
temperature  was  formulated  according  to  Holton  (1972)  as 


L  q 

0^  =  0  EXP  -^-^  ,  (4.9) 

sE  C   T 

P 


and  the  moist  convective  adjustment  was  designed  to  ensure 
that  the  moist  adiabatic  lapse  rate  was  not  exceeded 

(90   /3p  <  0)  . 

S  E 

The  heat  release,  Q  ,  as  a  result  of  moist  convective 

xm 

adjustment  was  included  through  knowledge  of  the  change  of 
equivalent  potential  temperature  during  the  adjustment. 

By  taking  the  natural  logarithm  of  equation  (4.9)  it 
is  found  that 


L  q 

In  0  „  =  In  0  +   °  *  (4.10) 

sE  C   T 


In  a  model  using  pressure  coordinates,  a  constant  pres- 
sure surface  may  be  assumed,  so  by  differentiating  both 
sides  of  equation  (4,10)  and  recalling  the  definition  of 
potential  temperature  gives 


d0  „    Jm    dq   L     L   q   dT 

sE    dT       s   c      c   s  ,  -  , , , 

-0— =T-  +  ^r^ ~2 — •  (4-11} 

sE  p  T 


Solve  the  Clausius  Clapeyron  equation  (4.6)  for  dq  , 

s 

substitute  in  equation  (4.11),  and  solve  for  dT.   The 
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solution  is  found  to  be 


dT  = 


d9  „ 
sE 

2 

6    il(1  +  IiLJL_   _  _^s 

SE   T      R  C  T2    V 
v  p 


)]  , 


(4.12) 


30 

where  dO    is  the  change  in  0    needed  to  adjust  = — 
sE  3       sE  J     dp 


sE 


=  0. 


Since  the  atmosphere  remains  saturated  throughout  this 
process  of  moist  convective  adjustment,  the  Clausius- 
Clapeyron  equation  was  again  used  to  obtain  the  correspond- 
ing dq  or  actually  dq  .   Once  again  the  changes  were  then 
applied  to  the  potential  temperature  and  specific  humidity 
fields  of  the  model  accounting  for  the  term  Q       of  equation 

(2.2)  and  M   of  equation  (2.3). 
c 

A  cumulative  amount  of  precipitation  was  then  found 
using  the  relation 


--*/ 


<5q  dp 


(4.13) 


in  agreement  with  the  changes  in  specific  humidity  due  to 
condensation  and  moist  convective  adjustment.  Precipita- 
tion was  expressed  in  inches  of  rainfall  in  the  model. 
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V.   NUMERICAL  SOLUTIONS 

All  numerical  results  to  be  shown  use  the  following 
values  for  the  constant: 


At  =  540  sec    p   =  100  cb    p   =  60  cb    H  =  80  cb 

*o  *m 

_2 
Ap  =  1.78  cb     Y  =  1800  km   Ay  =  60  km    g  =  9.81  m  sec 

f   =  10~4sec~    D  =  10~5sec-1  6  =  300°K 

o 

0  (pQ)  =  270,  280,  290,  300°K    C   =  1003  joules  kg~1°K~1 

R   =  461  joules  kg"  °K~     R  =  287  joules  kg"  °K~ 

6  -1 

K  =  R/C     L   =  2.5  x  10   joules  kg        T   =  273°K 
p     c  o 

a  =12.56°K       RH  =  50,70,85% 

3®I  -1      9®I  -1 

■= =  4°K  km    ,   (-r )   =  -.53448,  .6,  ,8°K  cb 

dz  dp   p 

m 

4     5   2-1 
A   =  A„  =A   =3x10,10   m  sec 
m     0     q 

2    -1 
C   =  0,  .001409,  .0001761  cb  sec     (which  corresponds  to 
m  e 

2    -1 
0,  10,  1.25  m  sec    in  height  coordinates) 

2    -1 
CQ  =  0,  .001409,  .0001761  cb  sec    (which  corresponds  to 
o 

2    -1 
0,  10,  1.25  m  sec    in  height  coordinates) 

2-1  2-1 

C   =0,  cb  sec     (which  corresponds  to  0  m  sec    in  height 

coordinates ) 

where  more  than  one  indicated  value  means  that  the  constant 

was  varied  in  different  experiments.   Whenever  possible, 

the  values  are  the  same  as  those  used  in  W73. 

In  this  section  numerical  f rontogenesis  solutions  with 

moisture  are  examined  at  t  =  4  days  and  compared  to  the 

corresponding  dry  atmosphere. 
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As  can  be  seen  from  Table  I,  26  experiments  were  per- 
formed on  five  cases,  cases  selected  from  Williams  (1973). 
Case  1  contains  only  horizontal  turbulent  diffusions  of 
momentum,  heat,  and  moisture  while  cases  2  through  5  con- 
tain both  horizontal  and  vertical  turbulent  diffusions 
of  these  parameters.   Values  chosen  for  these  turbulent 
diffusion  coefficients  were  computed  to  correspond  with 
values  used  in  Williams  (1973). 

In  Table  I,  0  (p  )  is  the  reference  potential  tempera- 
ture (°K)  as  introduced  in  equation  (3.7).   The  condensa- 
tion interval  indicates  the  number  of  time  steps  between 
successive  applications  of  condensation  and  moist  convec- 
tive  adjustment  to  the  model.   Turbulent  diffusion  coeffi- 
cients are  then  given  as  previously  defined  followed  by 
the  initial  relative  humidity  (percent)  and  the  initial 
static  stability (°K/cb) . 

Experiments  were  first  performed  with  a  reference 

potential  temperature  (9  (p  ))  of  300°K  to  directly  compare 

I   o 

with  the  dry  solutions  of  W73.   This  corresponds  to  a  tem- 
perature range  of  62-104°F  at  the  surface  along  a  North- 
South  axis.   Therefore,  experiments  were  also  conducted 
with  0  (p  )  =  270, 280, 290°K  to  simulate  physically  realis- 
tic situations.   Changes  to  the  dry  atmosphere  with  the 
different  references  were  almost  negligible,  but  due  to 
the  exponential  function  in  the  Clausius-Clapeyron  equa- 
tion important  differences  were  observed  in  the  moist 
atmosphere . 
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Experiments  were  also  conducted  using  different  values 
for  initial  relative  humidities  and  initial  static  stabili- 
ties o 

A  brief  tabulation  of  experimental  results  is  given  in 
Table  II,  comparing  widths  of  frontal  zone,  frontal  slopes, 
vertical  motions  and  precipitation. 

A  reasonable-  measure  of  the  width  of  the  frontal  zone 
is  given  by 


,  -  8  (-Y,s,o)  -  6  (Y,s,o) 

[Hi 

1 3y ' max 


(5.1) 


S8 
where  -^ —  is  approximated  by  a  one-sided  difference.   This 
dy 

width  (km)  is  indicated  at  three  levels  for  each  experiment 

in  the  table. 

The  North-South  location  of  |-tt— I     at  three  levels  is 

1 dy ' max 

y  (km)  and  is  used  as  a  position  indicator  of  the  frontal 
zone  where  |  y  ]  <^  1800  km. 

The  pressure  level  and  magnitude  of  maximum  vertical 
motion  is  then  given  followed  by  the  North-South  position 
of  maximum  precipitation  (km),  |y  |  £  1800  km,  and  the 
corresponding  amounts  of  cumulative  precipitation  over  the 
four  day  period. 

It  is  noted  from  Table  II  that  widths  of  frontal  zones 
between  corresponding  moist  and  dry  experiments  are  pri- 
marily effected  at  mid  to  upper  levels  with  only  small 
deviations  being  noted  at  the  surface.   The  expected  in- 
crease of  thermal  gradient  at  these  mid  to  upper  levels 
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due  to  latent  heating  in  the  moist  cases  is  seen  to  be 
the  cause. 

Comparisons  of  frontal  slopes  of  corresponding  experi- 
ments revealed  deviations  again  at  mid  levels.   The  slopes 
appear  more  vertical  throughout  these  levels  in  the  moist 
experiments.   Graphical  illustrations  will  later  be  pre- 
sented . 

Maxima  in  vertical  motion  seem  to  intensify  and  be 
found  at  higher  levels  along  the  frontal  zone  in  the  cor- 
responding moist  cases.   This  effect  can  also  be  directly 
tied  to  releases  of  latent  heating  in  this  region. 

Cumulative  precipitation  amounts  in  Table  II  are  ob- 
served to  be  primarily  functions  of  initial  relative 
humidity  and  reference  potential  temperature. 

Primary  emphasis  on  analyses  of  experimental  results 
was  on  contour  plots  from  the  IBM  360/67  computer.   These 
results  obtained  for  all  26  experiments  are  analyses  at 
t  =  4  days  of  the  fields  of  u,q,9'  and  w  where  G1  =  0  -  9 
and  9   is  the  horizontally  averaged  initial  potential  tem- 
perature.  The  entire  vertical  plane  is  shown  from  100  to 
20  cb  in  these  plots.   However,  the  outer  portion  of  the 
domain  between  |y|  =  1200  km  and  |y|  =  1800  km  is  not  shown, 
Also  note  that  isolines  of  specific  humidity  in  the  figures 
to  follow  have  different  intervals  dependent  upon  reference 
potential  temperatures  in  order  to  allow  for  proper  defini- 
tion of  the  field. 
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Due  to  the  large  number  of  contour  plots  obtained,  ex- 
periments 1,2,4,5,6,11,12,13,17,18  were  chosen  as  repre- 
sentative samples  for  detailed  discussion.   Experiments 
1,2,4,5,6  and  11  contain  only  horizontal  turbulent  diffu- 
sions of  momentum,  heat,  and  moisture  while  experiments 
12,13,17  and  18  contain  both  horizontal  and  vertical  turbu- 
lent diffusions  of  these  parameters. 

Figure  1  comparing  experiments  1  and  2  shows  the 
frontal  scale  at  99.1  cb  reaching  steady-state  at  approxi- 
mately t  =  2  days.   Little  deviation  is  seen  between  the 
wet  and  dry  experiment.   Figure  2  compares  the  frontal 
slopes  of  the  same  two  experiments.   As  mentioned  before 
slopes  are  seen  to  be  largely  the  same  with  the  only  signi- 
ficant deviation  being  at  midlevels  where  the  moist  slopes 
approaches  the  vertical  and  then  recurves  in  agreement 
with  the  dry  solution. 

Again  this  deviation  in  the  wet  experiment  is  probably 
due  to  the  release  of  latent  heating  at  these  levels. 
Similar  results  were  achieved  for  other  cases  with  more 
extreme  deviations  being  observed  at  higher  reference 
potential  temperatures. 

Contour  plots  for  experiments  1  and  2  are  shown  as 
figures  5  -  12.   Figures  5  and  6  show  distributions  of  the 
rotational  wind,  u  for  a  corresponding  dry  and  moist  experi- 
ment respectively.   It  is  noted  that  low  level  features  are 
quite  similar  in  the  easterly  flow.   However,  horizontal 
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and  vertical  wind  shear  is  more  intense  at  low  mid  to  mid 
levels  in  the  warm  air  in  the  moist  experiment.   It  is  note- 
worthy that  this  property  is  dominant  throughout  the  moist 
solution. 

Figures  7  and  8  contrast  values  of  9 '  throughout  the 
domain.   Noted  is  the  more  intense  gradient  of  potential 
temperature  through  the  mid  levels  and  the  pocket  of  warm 
air  at  mid  levels  most  likely  due  to  release  of  latent 
heating  in  the  moist  experiment. 

Fields  of  specific  humidity  are  contrasted  for  experi- 
ments 1  and  2  in  figures  9  and  10.   In  the  dry  experiment 
moisture  is  merely  advected  and  horizontally  diffused 
throughout  the  domain  showing  an  interesting  though  physi- 
cally unrealistic  situation.   The  processes  of  condensation 
and  moist  convective  adjustment  are  only  applied  in  the 
moist  experiments.   By  comparing  the  isolines  of  specific 
humidity,  it  is  seen  that  the  processes  of  condensation 
and  convective  adjustment  are  indeed  reducing  the  moisture 
content  of  the  air  and  that  the  gradient  of  specific 
humidity  is  stronger  along  the  frontal  zone  in  the  dry 
experiment  where  no  condensation  is  taking  place. 

Vertical  motions  are  investigated  for  the  experiments 
in  figures  11  and  12.   Upward  vertical  motions  are  stronger 
and  occur  at  higher  levels  closer  to  the  frontal  zone  in 
the  warm  air  mass  for  the  moist  experiment.   Downward  ver- 
tical motion  in  the  cold  air  mass  was  also  slightly  inten- 
sified with  its  maximum  also  at  a  higher  level. 
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Figures  13  -  28  compare  the  dry  and  moist  fields  of 
u,  q,  9*  and  w  at  reference  potential  temperatures  of 
290,  300  and  270°K.   Figures  13  -  16  show  moist  fields 
(8  (p  )  =  290°K).   The  horizontal  and  vertical  shear  of 
u  in  figure  13  is  seen  to  be  stronger  and  to  be  extending 
to  higher  levels  than  in  the  previous  moist  experiment. 
Stronger  thermal  gradients  are  seen  to  extend  higher  in 
the  atmosphere  in  figure  14.   Latent  heat  release  occurs 
at  higher  levels  and  to  greater  extent  as  is  seen  in  the 
upper  level  isentropic  pattern  to  the  warm  air  side  of 
the  frontal  zone. 

As  is  predicted  by  the  Clausius-Clapeyron  equation 
(3.18),  the  moisture  values  (figure  15)  are  significantly 
greater  than  in  the  previous  experiments.   Again  note 
the  decreased  gradient  of  moisture  along  the  frontal 
zone  in  the  moist  experiment. 

Vertical  motions, shown  in  figure  16,  are  also  somewhat 
intensified,  occur  at  higher  levels  and  have  slightly 
more  intense  gradients  along  the  discontinuity  in  the  warm 
air.   The  downward  vertical  motions  in  the  cold  air  are 
also  somewhat  increased. 

Figures  17  -  24  contrast  dry  and  moist  experiments 
with  a  reference  potential  temperature  of  300°K.   Little 
difference  is  seen  between  the  dry  experiment  and  the 
earlier  dry  fields  based  on  280°K  except  for  a  northward 
displacement  of  the  fields  at  the  colder  reference  poten- 
tial temperature.   Due  to  the  exponential  nature  of  the 
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moisture  parameter,  moist  fields  showed  significant  changes 
in  the  following  experiments. 

The  rotational  wind  u,  as  is  shown  in  figures  17  and 
18,  has  much  stronger  horizontal  and  vertical  shear  in  the 
warm  air  at  low  mid  to  upper  levels  than  in  previous  moist 
experiments  and  even  has  an  area  of  easterly  winds  at  upper 
levels  in  the  warm  air.   The  air  is  much  warmer  in  the 
upper  levels  due  to  latent  heat  release,  as  indicated  by 
figures  19  and  20.   Thermal  gradients  along  the  frontal 
zone  in  the  warm  air  are  also  stronger  in  the  mid  and  upper 
levels . 

Gradients  in  the  moisture  field  are  greatly  intensified 
and  specific  humidities  are  in  excess  of  28  g/kg  at  low 
levels  to  the  South.   Recall  that  the  "dry"  q  field  (figure 
21)  is  distributed  due  to  advective  and  horizontal  diffusive 
processes  only.   Perturbations  in  the  moist  experiment 
(figure  22)  are  probably  due  to  the  larger  quantities  of 
latent  heat  release  observed  in  this  experiment. 

Vertical  motions  (figures  23  and  24)  exceed  4  x  10 
mb/sec  in  the  warm  moist  air.   These  features  occur  at 
upper  levels  and  a  cell  containing  downward  motion  is 
observed  at  these  upper  levels  to  the  South  and  seems  to  be 
somewhat  due  to  the  impermeable  upper  boundary. 

Fields  with  a  reference  potential  temperature  of  270°K 
are  finally  examined  in  figures  25  -  28.   Due  to  minimal 
moisture  content  at  these  lower  temperatures  the  fields 
are  quite  comparable  to  dry  experiments  previously  shown. 
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Shear  in  the  u  field  (figure  25)  and  thermal  character- 
istics (figure  26)  are  also  quite  similar  to  dry  experi- 
ments.  Specific  humidity  (figure  27)  reaches  a  maximum 
of  5  g/kg  allowing  for  only  a  small  possible  latent  heat 
release.   Finally,  vertical  motions  (figure  28)  occur 
mainly  at  the  lower  levels  and  are  mere  closely  tied  to 
dry  experiments. 

Further  experiments  with  case  1  are  included  in 
Table  I.  Due  to  the  repetitive  nature  of  the  contour  plots 
obtained,  they  are  not  included  as  figures,  but  some  fur- 
ther discussion  is  appropriate  at  this  time. 

Experiments  3  and  8  were  run  with  varying  initial  rela- 
tive humidities.   Higher  initial  humidities  do  lead  to 
more  latent  heat  release  with  its  increased  thermal  gradi- 
ents, wind  shear,  and  vertical  motion  in  the  warm  air  but 
did  not  seem  to  warrant  further  experimentation. 

Experiments  9  and  10  were  run  with  varying  static  sta- 
bilities but  due  to  extensive  similarities  in  the  solu- 
tions, further  experiments  were  not  required. 

The  condensation  interval  was  generally  five  time  steps, 
The  iterative  procedure  of  returning  relative  humidities 
to  100%  worked  equally  well  at  everj  ten  time  steps,  never 
exceeding  two  iterations.   However,  five  time  steps  were 
finally  chosen  to  compromise  computer  time  and  physical 
reality.   Only  minor  changes  were  observed  in  the  contour 
plots  of  u,  q,  6'  and  w. 
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Experiments  that  follow  contain  both  horizontal  and 
vertical  turbulent  diffusions  of  momentum,  heat,  and 
moisture.   As  in  the  horizontal  case  previously  discussed, 
values  of  the  diffusion  coefficients  were  computed  to 
correspond  with  those  used  in  Williams  (1973). 

Figure  3  compares  the  frontal  slopes  for  the  corres- 
ponding dry  and  moist  experiments  12  and  13.  Once  again, 
slopes  are  almost  identical  with  the  exception  of  at  mid- 
levels,  where  in  the  moist  experiment  maximum  potential 
temperature  gradient  again  shifts  toward  the  warm  air  and 
then  recurves  to  agree  with  the  dry  experiment. 

Figure  4  for  the  same  experiments  shows  frontal  scale 
at  99.1  cb  versus  time  with  steady-state  being  obtained 
at  approximately , t  =  2  days.   Oscillations  arise  from  the 
development  of  the  surface  boundary  layer,  which  is  not 
present  in  the  initial  conditions  as  was  described  in 
Williams  (1973).  Very  little  difference  is  observed  between 
the  moist  and  dry  solutions  at  low  levels  which  reinforces 
the  idea  that  the  model  remains  essentially  unchanged  at 
these  levels  with  addition  of  moisture. 

Contour  plots  are  then  shown  as  figures  29  -  36  for 
experiments  12  and  13.   Compared  to  the  corresponding  dry 
case,  horizontal  and  vertical  wind  shear  is  again  increased 
especially  at  low  mid  to  upper  mid  levels  in  the  warm  air 
(figures  29  and  30) .   The  thermal  pattern  (figures  31  and 
32)  shows  increased  gradient  at  the  low  mid  to  upper  levels 
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in  the  warm  air  and  the  area  of  maximum  latent  heat  release 
is  seen  in  the  isentropic  pattern.   Vertical  influences 
for  the  moist  experiment  seem  somewhat  more  pronounced  for 
a  reference  potential  temperature  of  280°K  as  compared  with 
similar  experiments  containing  only  horizontal  diffusion. 

Figures  33  and  34  compare  the  fields  of  specific 
humidity.   The  changes  continue  to  be  a  reduction  in 
gradient  along  the  frontal  zone  in  the  prime  area  of  latent 
heat  release  in  the  moist  experiment.   Little  change  is 
observed  throughout  the  rest  of  the  domain. 

Vertical  motion  fields  (figures  35  and  36)  for  these 
same  experiments  show  low  level  maxima  of  1.5  x  10   mb/sec 
in  the  warm  air  but  an  additional  more  intense  maximum  of 
vertical  motion  is  found  at  mid  and  mid  upper  levels  along 
the  frontal  zone  in  the  moist  experiment. 

Case  3  includes  a  vertical  thermal  turbulent  diffusion 
coefficient  eight  times  as  large  as  in  case  2  previously 
discussed.   From  Williams  (1973)  a  decrease  in  frontal  tilt 
is  expected  along  with  an  increase  of  frontal  scale  due  to 
the  change  in  diffusion  coefficients.   This  can  be  observed 
by  comparing  appropriate  figures  of  the  dry  atmosphere. 
Other  changes  to  the  fields  will  most  likely  be  due  to  in- 
clusion of  moisture.   At  any  rate  our  comparisons  will 
remain  as  before,  a  comparison  between  the  corresponding 
dry  and  moist  experiments. 

Figures  37  and  38,  comparing  experiments  17  and  18, 
show  again  the  increase  of  horizontal  and  vertical  wind 
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shear  at  low  mid  to  upper  mid  levels  to  slightly  greater 
extent  than  in  the  previous  moist  experiment,  experiment 
30. 

The  isentropic  analyses  (figures  39  and  40)  show  the 
tightening  of  potential  temperature  gradient  at  low  mid 
to  upper  levels  once  again.   However,  it  is  now  observed 
that  the  increased  thermal  diffusion  has  smoothed  somewhat 
the  perturbation  in  the  isentropes  in  the  area  of  maximum 
latent  heat  release.   This  added  vertical  thermal  diffusion 
seems  to  aid  in  distributing  the  isolated  heating  induced 
"pocket"  that  had  been  noted  previously. 

Figures  41  and  42  compare  specific  humidity  fields  of 
the  same  experiments  showing  the  weakening  of  gradients  at 
low  mid  to  upper  levels  due  to  condensation. 

The  fields  of  vertical  motion,  figures  43  and  44,  are 
similar  in  results  to  the  previously  discussed  experiment 
containing  horizontal  and  vertical  turbulent  diffusions 
except  for  some  smoothing  of  features  due  to  the  added 
vertical  thermal  diffusion  of  the  present  experiment.   Once 
again  the  maxima  noted  in  the  dry  experiment  18  is  main- 
tained with  an  added  maximum  in  vertical  motion  in  the  mid 
level  warm  air  along  the  frontal  zone  in  the  moist  experi- 
ment . 

After  an  analysis  of  the  before  mentioned  experiments, 
it  becomes  obvious  that  an  observational  pattern  has 
developed  and  that  further  discussion  using  the  line  of 
attack  would  be  largely  repetitive. 
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Careful  analysis  of  the  contour  plots  of  the  remain- 
ing of  the  26  experiments  contained  essentially  the  same 
features  as  those  discussed  herein. 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  numerical  model  of  Williams  (1973)  is  modified 
in  this  study  to  include  moisture.   The  numerical  solutions 
become  quasi-steady  within  1-2  days  as  in  W7 3 .   Intensifi- 
cation of  the  baroclinic  zone  at  mid  to  upper  levels  seems 
to  be  the  significant  addition  to  the  model  due  to  the  in- 
clusion of  moisture.   This  intensification  of  the  baro- 
clinic zone  at  the  mid  to  upper  levels  is  consistently 
manifested  by  the  increases  of  thermal  gradient,  wind 
shear,  and  vertical  motions  in  the  warm  air  along  the 
frontal  zone  which  are  due  to  latent  heat  release.   Second- 
ly, it  should  be  noted  that  the  model  is  essentially  un- 
changed at  low  levels  with  the  inclusion  of  moisture. 
Thirdly,  the  model  with  moisture  shows  dependence  upon 
the  reference  potential  temperature,  6  (p  )  chosen  due 
to  the  exponential  relationship  between  moisture  and  tem- 
perature as  is  indicated  by  the  Clausius-Clapeyron  equa- 
tion (3.18).   It  is  recalled  that  the  dry  model  is 
relatively  insensitive  to  choice  of  this  reference  tem- 
perature.  The  moist  model  illustrates  this  exponential 
thermal  dependence  in  the  various  experiments  at  the 
reference  potential  temperatures  of  270,  280,  290,  and 
300°K. 

In  future  studies  using  Williams  model,  attention 
should  be  directed  towards  transformation  to  In  p 
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coordinates  for  a  better  depiction  of  upper  level  gradi- 
ents.  Upper  level  boundaries  might  also  be  changed  and 
raised  to  include  a  stratosphere  and  finally  a  more 
sophisticated  convective  scheme  might  be  applied  to  the 
model  such  as  that  of  Arakawa  and  Schubert  (1973). 
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Figure  1.   Frontal  scale.   Experiments  1  and  2, 
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Figure  2.   Frontal  slope.   Experiments  1  and  2. 
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Figure  3.   Frontal  scale.   Experiments  12  and  13. 
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Figure  4.   Frontal  slope.   Experiments  12  and  13. 
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Figure   14.      6'(°K).      Exp.    4.  Y(km) 
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Figure    16.    w(xl0     mb/sec) .      Exp. 4. 
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Figure    18.      u(m/sec).      Exp.    6. 


0 
Y(km) 


■300  -600  -900  -1200 


58 


20- 

2  K*rinT*'         " 

» 

j 
40- 

?   60- 

a. 

\Wl 

80- 

xVxxvA     / 

• 

\^\^^v  ( 

100- 

^^■^itlkV 

12*00  900  600  300  °  ~300  -600 

Figure   19.      9'(°K).      Exp.    5.  Y(km) 


-900     -12 


00 


1200     900      600      300      0 
Figure  20.   e'(°K).   Exp.  6.       Y(km) 


•300     -600 


i       r 

900    -1200 


59 


20- 


40- 


v    60- 


80- 


100- 


12'00  900  600  3&0  0 

Figure  21.      q(gAg).      Exp.    5.  Y(km) 


-300 


-600  -900  -1200 


1200  900  600  300  0  -300  -600 

Figure   22.      q(g/kg).      Exp.    6.  Y(M 


-900  -1200 


60 


20- 


40- 


£   60- 


80- 


100- 


'  12*00 


900  600  3&0  0  -300 

Y(lcm) 

Figure   23.    w(xl0     mb/sec) .    Exp.    5. 
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Figure   24.    w(xl0     mb/sec).      Exp.    6. 
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Figure    28.    w(xlO      mb/sec) .      Exp.    11. 
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Figure   30.      u(m/sec).      Exp.    13.        Y(km) 
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Figure   32.      e'(°K).      Exp.    13.  Y(km) 
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Figure    34.      q(g/kg).      Exp.    13.  Y(km) 
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Figure   36.    w(xl0     mb/sec) .      Exp.    13. 
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Figure   38.      u(m/sec).      Exp.    18 
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Figure   39.      Q' (°K) .      Exp.    17.  Y(fcm) 
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